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ABSTRACT 



We analytically work out the long-term variations caused on the motion of 
a planet orbiting a star by a very distant, pointlike massive object X. Apart 
from the semi-major axis a, all the other Keplerian osculating orbital elements 
experience long-term variations which are complicated functions of the orbital 
configurations of both the planet itself and of X. We infer constraints on the 
minimum distance dx at which X may exist by comparing our prediction of the 
long-term variation of the longitude of the perihelion cjj to the latest empirical 
determinations of the corrections Aw to the standard Newtonian/Einsteinian 
secular precessions of several solar system planets recently estimated by indepen- 
dent teams of astronomers. We obtain the following approximate lower bounds 
on dx for the assumed masses of X quoted in brackets: 150 — 200 au (Mars), 
250 - 450 au (0.7 m®), 3500 - 4500 au (4 mj^p). 

Subject headings: Planets • Planet X • Solar System perturbations 
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Introduction 



Does not any other planet of the Sun exist in addition to those already known? Have 
not the remo te peripheries of the solar system saved any furthe r big surprises for us? It is 
a long tim e (iLowelllllQlSl: |Pickerindll928Ul931t ISchuettelll949[) that such questions— still 
well alive (ILykawka and Mukail l2008l : llorid l2010l : iMatese and Whitmird l201ll : [Fernandez 
20 111 )— are well rooted in the astronomical community, having often resurface d in di fferent 
contexts and with changed forms;. For a recent, popular review, see ISchillinig ( l2009l ). Here 
we limit ourselves to recall that 



The hypothesis of a stellar-like companion object (Nemesis) orbiting the Sun w as 
postu lated for explaining terrestrial extinction per iodicitv (IRaup and Sepkoskil 
19841). thought to be mediated by cornet showers (IWhitmire and JacksonI Il984j: 
Davis et all 1 1984 lHillslll984l : lHut[ll984l : IVandervoort and Sathei]ll993l : lMullei]l2002h . 
A recent critical review about the evidence for and agai nst astronomi c al irn pacts 
on climate change and mas s extinctions can be found in iBailer- JonesI fl2009h : see 
Melott and Bambachl (120101 ) for recent developments of the investigations on the role 
of Nemesis in such phenomena. 

Anomalies in the distribution of the aphelia and the orbital elements of the comets 
resident in the outer region of the Port cloud led to postulating the existence of a 



Sun-bound Jovian mass body flMatese et al.lll999l ) 



among a putative Planet X and comets are, e.g. 



Matese et al. ( 


1986 


); 


Murrav 



a Jupiter-type object, now named Tyche (IMatese and Whitmird 1201 ll ) and not to 



be confused with Nernesis si nce it would not be able to induce comet s torms, see 



torn . 



Matese and Whitmird (120111 ). Another recent analysis can be found in iFernandez 



postulated ( 


Vlaran et al. 


1997; 


CoUander-Brown et al. 


2000; 


Brunini and Melita 


2002 


Melita et al. 


2003. 


2004: 


Matese et al. 2005: Gomes et al. 


2006: 


Lvkawka and Mukai 


2008 


), with changeable fortune, to explain several features pertaining the dynamical 



history and the architecture of the Trans-Neptunian Belt. See iLykawka and Mukai 



(12008! ) for recent developments of such a scenario. 



A putative Planet X was also invoked 


Gunn 


1970: Rawlins 


1970: Seidelmann 


1971: ^ 


v^an Flandern 


1981; 


Anderson and van Flandern 1982: 


Anderson 1987: 


Seidelmann and Harrington 1988: Harrington 


L988: Gomes and Ferraz-mello 1988: 


Gomes 


1989: Powell 


1989; 


van Flandern 


1991) to explain certain seeming irregularit 



i n the orbital motions of Uranus and Neptune w hich show e d up in ancient data 



(jBruninilll992l ): such an issue was later settled by IStandishI (119931 ) with a re-analysis 



of the observations by including more recent radiometric points from the Voyager 2 
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automatic spa cecraft w h ich al lowed for a more accurate determination of the mass of 
Neptune. See IStandishI ( 1l996l ) for issues related with the orbital residuals of Pluto, 
likely of non-dynamical origin. 



HarrisonI (119771 ) suggested that, as an explanation of the peculiar properties of certain 
pulsars with anomalously small period derivatives, the barycenter of the solar system 
is accelerated, possibly because of a hithe r to un detec ted companion star of the Sun i n 
a bound or open orbit. See also ICowlingl ( 119831 ) and IZakamska and Tremaind (120051 ) 
for further studies on such a topic. 



loriol (120101 ) found that a distant, pointlike massive object in the outer regions of the 
solar system may explain the ar iomalous retr ograde precessioiJl] of the perihelion of 
Saturn preliminarily estimated ( jPitjeval 120081 ) by analyzing some radiotechnical data 
points from the Cassini spacecraft exploring the Saturnian system. 



Here by means of the expression "Planet X" we broadly denote a pointlike object, having a 
mass approximately as large as that of Mars or larger, up to typical values of brown dwarfs 
{nix ~ 80 mjup) or even red dwarfs (mx ~ 0.5 Mq) according to certain scenarios, located 
at large distances from the Sun, and for which no direct observational evidence, based on 
the detection of electromagnetic radiation of different wavelengths emitted or reflected by 
it in natural processes, is (yet?) available. 

We also note that investigations on P lanet X have di r ect connections with funda mental 
physics as well since it has been shown ( iMilgronj l2009l : iBlanchet and NovakI 1201 ll ) 
that the action of a distant body located towards the Galactic center is dynamically 
equivalent to that of the External Field Effect (EFE) in the planetary regions of the solar 
system within the f r amew ork of the Modified Newtonian Dynamics (MOND). Moreover, 
Foot and Silagadzd (120011 ) considered the possible existence of solar system's planets made 
of non-annihilating mirror matter, which is one of the candidates for the non-baryonic Dark 
Matter which is believed to constitute about 22% of the matter content of the Universe 



( iBertone et al.l 120051). For o t her r esearches on the existence of mirror matter in the solar 
system, see iFoot and Mitral ( l2003l ): a general overv iew on the mirror matter scenario and 



its observat ional implications can b e foun d in. e.g.. iFootI (20021). For earher studie s on such 



a topic, see iBlinnikov and Khlopovl ( 119831 ): iKhlopov et al. 



(Il99lf ): iKhlopovl ( 119991 ). Finally 



the subtle perturbing effects due to a putative X may be important in realistically assessing 
the error budget in several high-precision tests of standar d post-Newtonian gravity proposed 
or to be performed in the solar system arena (jloridl201ll ). 



In this paper, we use recent and accurate observational determinations of the motion 
of some planets of the solar system to put tighter constraints on the location of Planet X 



^Anyway, more recent data analyses (|Pitjeva]l2010l ) yield values for such an anomaly which 
are statistically compatible with zero. See Section [31 
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with respect to those existing in hterature. The plan of the paper is as follows. In Section 
|2]we analytically work out the effects that a distant body X would cause on the orbital 
motion of a closer and (relatively) fast-moving planet. More specifically, we adopt the 
Lagrange perturbative scheme to calculate, without simplifying approximations concerning 
the orbital configurations of both the planet and of X, the long-term variations of all the 
Keplerian orbital elements of the perturbed planet. In Section [3] we turn to the latest 
observational determinations of the planetary motions. In particular, we compare our 
theoretical prediction for the long-term variation of the perihelion with the most recent 
determinations of the corrections to the standard precessions of the perihelia of the major 
bodies of the solar system obtained by independent teams of astronomers with the latest 
ephemerides. Our goal is to look at the minimum distance at which Planet X may exist as 
a function of its position in the sky. In Section H] we review the past and future sky surveys 
and compare our bounds on the minimum distance of X to their obtained and expected 
results. Section \5\ is devoted to the conclusions. 



2. Analytical calculation 



Let us consider a remote, point-like object X of mass mx and located at distance rx 
from a central body of mass M. Its action on a closer test parti cle orbiting M at distance r 
can be modelled as due to the following quadrupolar potential (iHogg et al.lll99ll ). accurate 
up to terms of order 0{r'^/r'^), 



3 r-l 



where 



Gm 



X 



X 



(2) 



'X 



is the tidal parameter of X, while 1 = {Ixyl'yJz} is a unit vector directed towards 
X determining its position in the sky. We purposely will not adopt any specific 
parameterization for it in order to make our results as general as possible. In eq. ([T]) 
r = {x, y, z} refers to the perturbed test particle. 

In order to compute the effects of eq. ([1]) on the orbital motion of a test particle 
we shall use standard pert urbative techniques. The use the eccentric anomaly E 
( iMurray and Dermottlll999l ) as fast variable of integration turns out to be more convenient 
to make the forthcoming integration more tractable. The average of eq. ([T]) over one orbital 
revolution of the perturbed test particle is 



n 



Ux dt 



32 



U{I,n,uj;l 



(3) 



- 6 - 



with 
U -- 



- (2 + 3e2) (-8 + 9/^ 



6/^) — 120e^ sin 2a; {l^ cos f2 + sin [1^ sin /- 



+ cos / {ly cos - sin Vi)] - ISe^ cos 2a; [3 (/^ - /2) cos 2VL + 2 [ll + - 2lfj sin^ /- 
- 4/^ sin 2/ {ly cos 1^ - sin fi) + Qljy sin 2^]] - 6 (2 + Se^) [ (/^ _ /2) cos 2Vt sin^ J+ 
+ 2U sin 2/ (/y cos Vt - sin 1]) + 2/^/^^ sin^ / sin 2n] - 3 cos 2/ { (2 + 36^) [ll + - 2lfj + 
+ 5e2 cos 2uj [ [ll - ll) cos 21] + 2/^/j, sin 2Vl\ ] . 

(4) 

In eq. ([3]) and eq. (jlj) a, e, /, fi, u are the semi-major axis, the eccentricity, the inchnation of 
the test particle's orbit to the reference {X, Y} plane, the longitude of the ascending node 
and the argument of pericenter, respectively; n = ^/GMjo?' is the unperturbed Keplerian 
mean motion connected to the unperturbed orbital period Pb by = 27r/Pb- Note that eq. 
([3]) and eq. (HI) are exact: neither approximations in e nor in J were used. In the integration 
/ was kept fixed over one orbital revolution of the perturbed test particle, as it is reasonable 
given the assumed large distance of X with respect to it. 



The standard Lagrange planetary equations (IBertotti et al.ll2003l ). in which the 
disturbing function, i.e. (f/x) in our case, appears, can be used to straightforwardly derive 
the long-term variations of all the standard six Keplerian osculating orbital elements of the 
perturbed particle. We also considered the longitude of the pericenter w = ^ -\- u since its 
secular precession is one of the parameters usually estimated by the astronomers wh en they 
fit the dynamical force models of their ephemerid es to huge planetary data records (jPitjeva 
20051 . 12OO8I . I2OIOI : iFienga et al.ll2010l : lFiengall2010h . 



From the Lagrange planetary equation for the semi-major axis (IBertotti et al.ll2003l ). 
it can be noted that a does not change since eq. (jl]) does not contain the mean anomaly 
M.. Instead, all the other Keplerian orbital elements experience non- vanishing long-term 
variations. They are 



with 
Z = 



de 15/CxeVl — 



dt 



IQn 



S i 1,^,00; I 



(5) 



— 8/^ cos 2u sin I {1^ cos Q + ly sin Q) + 4 cos / cos 2uj [—2lxly cos 2Q+ 

{ll - ll) sm2n] + sin 2a; [{ll - Py) (3 + cos 2/) cos2fi + 2 (/^ + - 21^) sin^ I- 

41^ sin 21 {ly cos Q — l^ sin ^2) + 2l,j.ly (3 + cos 21) sin 2VL\ . 

(6) 
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dl 



3/C, 



Z [1,^1, CO] l) , 



(7) 



dt 4 a/1 - e^n 
in which 

X = [/^ cos / + sin J (— cos r2 + Z^. sin fi)] {5e^/^ sin Jsin 2cj+ 

+ 5e^ cos / sin 2u {ly cos f2 — sin f2) + [2 + (3 + 5 cos 20;)] {l^ cos VL + ly sin fi) } 



with 



with 
V 



dn 

'dt 



A/" = 3 CSC / [/^ cos J + sin J (— Zj^ cos + /j. sin fi)] {— 2/^ sin /+ 
+ cos / [—2 + (—3 + 5 cos 2uj)] {ly cos — Ix sin f2) + 
+ [/^ (—3 + 5 cos 2uj) sin J — 5 sin 2u) {l^ cos 1] + /j^ sin f2)] } 



(9) 



(10) 



do; 

lit 



3/C 



X 



lenvT 



(11) 



= - (1 - e2) (-8 + 9/2 + 9/2 + 6/2) + (/2 + /2 _ 2/2) [_2 _ Se^ - 5 cos 21- 

+ 5 cos 2a; (-1 + 2e'^ + cos 2J)] + cos 2VL { [ll - Py) [-1 + 6e'^+ 

+ 5 (26^ - 3) cos 2u + 10 cos 21 sin^ u] - 20 (2 - e^) /a^/^^ cos / sin 2a;} + 

+ 21, cos Q {2ly cot / [(1 - e^) (-3 + 5 cos 2uj) + 10 cos 21 sin^ - 

- 5lx [2 — 3e^ — (2 — e^) cos2J] esc / sin 2a;} — /^ esc / sinf2 {Al^ cos / [(1 — 

- e^) (-3 + 5 cos 2a;) + 10 cos 21 sin^ a;] + 10/^ sin 2a; [2 - 36^- 

- (2 - e^) cos 2/]} + 2 sin 2Q [-5 (3 - 2e'^) IJy cos 2a; + IJy (-1 + Qe^+ 
+ 10 cos 2/ sin^ a;) + 5 (2 - e^) (/2 - /2) cos J sin 2a;] . 



X 



dw /C 
cit ~ 128nVl - e2 



^ /,fi,a;;n +n(I,n,uj;l 



(12) 



(13) 



where 

-Wt^) = 8-9/^-9/2_6/2_io(/2 + /2_2/2),og2a;sin2/+ 

+ 8lz cos Q sin / [ly cos / (—3 + 5 cos 2u) — 5lx sin 2uj] + 

+ cos 2Q [-3 {ll - ll) (5 cos 2u + 2 sin^ /) - 401 Jy cos / sin 2u] - 

— 81 z sin / [/^. cos / (—3 + 5 cos 2a;) + 5ly sin 2uj] sin f2+ 

+ 2 [-3/^/j^ (5 cos 2cj + 2 sin^ /) + 10 {ll - ll) cos / sin 2uj] sin 2Vt- 

- cos 21 {3 + /2 _ 2/2) + 5 cos 2a; [{ll - ll) cos 21] + 21 Jy sin 2fi] } , 

(14) 

and 

% = — 96 [/z cos / + sin/ (— Zj, cos f2 + /a; sin n)] {— 2Z2 sin /+ 

+ cos/ [-2 + (-3 + 5cos2a;)] (/j^cosfi - /x-sinfi) + (15) 
+ [Iz (—3 + 5 cos 2a;) sin / — 5 sin 2a; {l^ cos fl + ly sin f2)]} tan (//2) . 

with 

A = (7 + 3e2) (-8 + 9/2 + 9/2 + 6/2) + 120(1 + e2) sin 2a; (/^ cos f]+ 

+ /j^ sin Q) [Iz sin / + cos / {ly cos Q — l^ sin (])] + 15 (1 + e2) cos 2a; [3 {ll— 

- /2) cos 2fi + 2 (/2 + /2 - 2/2) sin^ / - 4/^ sin 2/ (/^^ cos n-k sin fi) + 
+ 61 Jy sin 2^] + 6 (7 + 3e2) [ (/2 - /2) cos 2^ sin^ / + 2/^ sin 2/ (/^ cos 

- Ix sin n) + 2/^/y sin^ / sin 20] + 3 cos 2/ { (7 + 3e2) {ll + 11^ 2lfj + 

+ 5(1 + e2) cos 2a; [{ll - ll) cos 2VL + 2/^/j, sin 2fi] } . 

It is not possible to simphfy such formulas by a — priori orienting one of the reference axes 
along I because the position of X is not known. In the case of the solar system, the reference 
{X, y} plane is typically the mean ecliptic at the epoch J2000.0, with the X axis directed 
towards the mean equinox at the same epoch. Our calculation holds for any orientation of 
the reference frame. 



(17) 
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Other calculations, performed in the fram ework of the lunisolar perturbations on 
artificial planetary satellites, exist in literature ( lKozailll959l : ICooklll962l : iGurfil et al.ll2007l ). 
They differ from ours in the choice of the orbital elements investigated, the calculational 
techniques adopted, and the representation of the unit vector of t he disturbing body. A 
straightforward comparison with our results is, thus, not possible. iKozail (jl959i ). by using 
a specific parameterization for the geocentric lunar unit vector / in terms of some orbital 
elements of it, adopted the Lagrange planetary equations for a,e, I,Q,u. Anyway, he did 
not average the disturbing function over one orbital revolution of the sat ellite, retainin g 
only the terms in which its mean longitude does not appear. Moreover, iKozail fll959h 
did not release explicit expressions for the variations of e, I at all, while those for Q and 
u are pre sumab l y inco mplete and difficult to evaluate since it is unclear how they were 
obtained. Cook ( 1962 ) used the more cumbersome Gaus^ equations for the variations of 
the osculating orbital elements to compute the long-term variations of a, e, / , fl, u + ^ cos / 
by means of the true anomaly / as fast variable of integration. Also ICookl ( 1l962l ) used 
a particul ar representa t ion o f the geocentric I, expressed in terms of its Keplerian orbital 



elements. iGurfil et al.l ( 120071 ). working in the framework of a perturbed Martian satellite, 
employed the satellite's mean anomaly with the Lagrange planetary equations in the 
integration over one orbital revolution for a, e, /, fl, io. They also adopted two specific 
parameterizations for I related to different choices for the reference {X, Y} plane. 

Our results can be used not only in view of the fact that estimated corrections Aw 
of the standard Newtonian/Einsteinian precessions of t he longitudes of the perihelia are 



nowadays ay ailable for all the eight planets and Pluto (lPitjevall2010t iFienga et al.l 12010 



Fiengall2010[ ). but also because it se ems that also the nodes will receive their due attention by 
the astronomers in the near future ( lFiengall2010l ). The first determination of the corrections 



to the ir standard secular precession for the first six planets can be found in iFienga et al. 
(l201l[ ): we propose to use them in further analyses. Moreover, eq. ([5])-eq. ( fT6l) can also be 



useful in severaH exoplanetary sc enarios. Indeed, in m any cases the discovered extrasolar 
planets have distant companions (IHoward et al.ll2010l ). or their existence is postulate d to 
explain certain observed features of the orbital motions ( Collier Cameron et al. 2010l ). 



3. Confrontation with the observations 

In Table [1] we quote the latest determinations of the corrections Ava to the standard 
Newtonian/Einsteinian secular precessions of the longitudes of the perihelia of all the major 
bodies of the solar system estimated by various authors with some of the most recent 
ephemerides. The effect of Planet X was not explicitly included in the dynamical force 



^He dubbed them as Lagrange's planetary equations. 
^See http://exoplanet.eu/ on the WEB. 
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Table 1: Estimated corrections Aw, in milliarcseconds per century (mas cty of the 
standard Newtonian/Einsteinian secular precessions of the l ongitudes of t he perihelia w 
of the eight planets and Pluto determined with the EPM 2008 f|Pitievall201oh . the INPOP08 
( Fienga et al.l 2010), and the INPOPlOa ( Fienga 2010) ephemerides. Concerning the val- 
ues quoted in the third c olumn from the left , they are those corresponding to the smallest 
uncertainties reported in (IFienga et al.ll2010l ). Note the small uncertainty in the correction 
to the precession o f the terrestrial perihelion, obtained by processing Jupiter VLBI data 
fiFienga et al.ll2010h . 



Planet 


Aw fPitieva 20101 


Aw fFienffa et al. 20101 


Aw fFienea 20101 


Mercury 


-4 ±5 


-10 ±30 


0.2 ±3 


Venus 


24 ±33 


-4 ±6 




Earth 


6±7 


0± 0.016 




Mars 


-7 ±7 


0±0.2 




Jupiter 


67 ±93 


142 ± 156 




Saturn 


-10± 15 


-10 ±8 


0±2 


Uranus 


-3890 ± 3900 


± 20000 




Neptune 


-4440 ± 5400 


± 20000 




Pluto 


2840 ± 4510 







models fitted to the observations, so that, in principle, the extra-rates of the perihelia in 
Table [U account for it. Here we shall compare our theoretical prediction of eq. ( 1T3|) to the 
corrections of Table [1] in order to infer upper bounds on /Cx? so to constrain the minimum 
distance dx at which X may be located for different values of its mass mx. 

By expressing the unit vector Z of X in terms of its ecliptic latitude /3x and longitude 
Ax and by equating w of eq. ( IT^ to Aw of Table [T] for each planet, it is possible 
to plot the minimum distance dy^ of X as a function of /3x, Ax for given values of its 
putative mass mx. In Figured] we plot such bounding surfaces obt ained from the accu rate 



extra-precess ions of the perihelia of the Earth, Mars and Saturn (IFienga et al 



2010 



Fienga! I2OIOI ) for mx = 4 mjup, correspondi ng to Tyche fjMatese and Whitmird 
mx = 0.7 m®, considered in the scenario by iLykawka and Mukail (120081 ) . and mx = mMars- 



2OIID 



It turns out that the bounds on (ix at high ecliptic latitudes are rather independent of 
Ax. For each of the perihelia there are a few sky locations, mainly close to the ecliptic, in 
which the dynamically inferred minimum distances of X are unrealistically small. Anyway, 
such positions differ from each other in such a way that a compensation occurs. Thus, 
a Tyche-like body cannot be located at less than about 3500 — 4500 au, an Earth-sized 
mass should not exist at less than about 250 — 450 au, while a rocky planet as big as 
Mars cannot be closer than approximately 150 — 200 au. The bounding surfaces by the 
Saturnian perihelion are generally smoother, with a milder excursion among the extremum 



Fig. 1. — Constraints on the minimum distance dx, in au, for rrix = 4 mjup (first row from 
the top), nix = 0.7 (second row from the top) and mx = "'TiMars (third row from the 
top) as a function of the echptic latitude Px and longitude Ax of X. The perihelia of the 
Earth (left column), Mars (middle column) and Saturn (right colu mn) were used acco rding 
to Ari70 = ± 0.016 mas cty~^ and Ap Mars = ± 0.2 mas cty~^ ( iFienga et al.ll2010l ). and 
Aiijsat = ± 2 mas cty"^ flFiengall201oh . 



values. They are rather close to those determined by the perihelion of the Earth, whose 
extra-precession is 2 orders of magnitude more accurate than for Saturn (cfr. Tabled]). 
Future improvements in the accuracy in determining the orbit of the ringed planet from 
continuous tracking of Cassini will straightforwardly yield more stringent constraints. It is 
possible to infer dx also for values of mx different from those used in Figure [H dubbed mx, 
by scaling the corresponding lower bounds dx by the multiplicative factor ^x = (i^x/ttt^xY^^- 

Our bounds are tighter than those obtained by_ llorid (l2009h with the correction to the 
precession of the perihelion of Mars estimated by iPitjeva (12005h . 
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Possibilities of direct detections 



Concerning the ability of several performed, ongoing or planned wide area sky surveys 
to directly reveal a putative Planet X from its emitted or reflected electromagnetic radiation 
of different wavelengths, let us recall the following. 

Limiting ourselves to optical wavelengths only, bodies as large as Jupiter or Neptune 
would have apparent visual magnitudes of about 20 and 23, res pectively, at 1000 au, so that 
the largest telescopes may image them like just very dim dots (jFernandezll201lt ). At 10000 
au it would impossible to detect them with any optical telescope available today (IFernandez 
201ll ). Actually, our tightest constraints from the Earth's perihelion pose a jovian-sized 
object at not less than approximately 3000 — 4000 au; the minimum allowed distance for 
a Neptune-like mass would be of no more than about 1110 au. Since low-temperature 
bodies emit more efficiently in the infrared, surveys operating at such wavelengths are more 
suitable to look for Planet X, especially those based in space because of the atmospheric 
extinction in the infrared. 

All the wide area sky surveys performed so far mostly explor ed regions close to the 
ecliptic. The ecliptic, all-sky optical survey by ( iTombaughl Il96ll ) concluded that the 
minimum distance at which an Earth-sized body coul d be located is 81 au; similar findings 
were reached by the nearly ecliptic optical survey by iKowall (119891 ). Our dynamically 
inferred lower bounds are about 3 — 6 times more stringent. The Spacewatch optical survey 



( iLarsen et al.l 120071 ). conducted within 10 deg from the ecliptic, was sensitive to Mars-sized 
objects out to 300 au and Jupiter-sized planets out to 1200 au; its negative findings, 
especially with respect to the detection of a Jupiter-like body, are well explained by our 
analysis since its minimum distance dynamically inferred from th e planetary perihelia is 
justmore than 1000 au for those ecliptic latitudes. According to iMatese and Whitmire 
( 120111 ). Planet X would have to be greater than 7—10 mjup and closer than 6000 — 25000 
au for a p ossible detection in th e infrared, ground-based Two Micron All Sky Survey 
(2MASS) ( ISkrutskie et al.l 120061 ). Our minimum distance for a body with mx = 7 mjup 
may be of the order of 1300 au in a few locations near the ecliptic, and 4800 au for high 
latitudes. iMatese and Whitmird ( 1201 ll ) note also that the negative searc h results of the 
InfraRed Astronomical Satellite (IRAS) surve }S (iNeugebauer et al.l 1 1984 ) suggest that an 
object with mx = 2 — 5 mjup must have a current minimum distance dx = 2000 — 10000 
au, respectively; our dynamically inferred bounds for the minimum allowed distance of such 
a X are in the range ~ 3000 — 4500 au for most of the positions in the sky. The all-sky 
synoptic Tycho-2 optical survey excluded the prese nce of a ra ain-sequence star above the 
hydrogen-burning limit within 1 pc= 2.06 x 10 au dHdlSQOOl). 



Certain rumors appe ared in 1983 (IO'Toolelll983l) in mass-media about IRAS and Planet 



X, indirectly caused by iHouck et al. 



spread by mass- media, appeared later ( iHouck et al.lll985l : ISoifer et al.lll987l ). 



19841). Clarificat i ons in specialistic literature, not 
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Moving to ongoing or forthcoming projects, the all-sky synoptic surv ey by the 
Panoramic Survey Telescope and Rapid Response System (Pan-STARRS) f JewittI 2003 ) 
would be able to reveal massive planets such as Neptune not beyond about 800 au, while 
a body with mx = 0.1 Mq would be undetectable for rx > 2000 au. Conversely, our 
dynamically inferred lower bounds for such kind of bodies in the ecliptic are up to ~ 1100 au 
and ~ 13400 au, respectively. A Ju piter-sized planet at 200 au over the whole sky could, 
in principle, be detected by GAIA ( iGaudi and Bloomll2005l ) with astrometric microlensing; 
actually, our lower limits for the admissible distance for such a body are globally more 
stringent. Also the putative existence of a X with mx = 3 mjup at 20000 au could be 
put on the test (iGaudi and Bloomll2005l ). Possible planets in the oute r regions ('>1000 
au) of the solar system may be revealed by the mesolensing technique (jPi Stefandl2008l ). 
Another future all-sky survey which may be useful in de tecting distant obje ct s far from the 



eclipti c is the Large Synoptic Survey Telescope (LSST) f lStubbs et al.ll2004l ). iBabich et al. 
( 120071 ) proposed to use observations of the Cosmic Microwave Background (CMB) to place 
constraints on the mass , distance, and size d i stribu tion of small objects in the Kuiper Belt 
and inner Oort Cloud. iMaris and Buriganal (120091 ) envisaged the possibility of using the 
data from the ongoing Planck mission to detect the thermal emission of pointlike objects 
of the solar system . The recently launched Wide-field Infrared Survey Explorer (WISE) 
( IWright et al.l 120101 ) should be able to detect a Neptune-sized object out to 700 au, which 
is, in fact, a distance smaller than our dynamically inferred bounds on dx for such a body. 
Instead, a Jupiter-mass object would be detectable out to 1 light year= 63239.7 au, where 
it would still be within the Sun's zone of gravitational control. A lar ger object of 23 Ju piter 
masses would be visible at a distance of up to 7 to 10 light years (ILakdawallal 120091 ) . 
i.e. 4.4 — 6.3 x 10^ au. Actually, WISE has recently discovered its first ultra-cold brow n 
dwarf, named WISEPC J045853.90+643451.9, at about 6 - 10 pc (iMainzer et al.l l20 111 ) 
corresponding to 1.2 — 2.0 x 10^ au. 



5. Summary and conclusions 

We analytically computed the long-term variations of all the Keplerian osculating 
orbital elements of a test particle orbiting a primary and acted upon by a pointlike, distant 
object X. We assumed that it is located at a much larger distance rx than the perturbed 
particle, so that the spatial position of X was kept fixed in the integration over one orbital 
revolution of the particle. With a perturbative, first-order calculation accurate up to terms 
of order 0{r'^ /r"^ we found that all the Keplerian osculating orbital elements, apart from 
its semi-major axis a, undergo long-term variations depending on the orbital geometries of 
both the perturbed and the perturbing objects in a complicated way. 

We compared our analytical prediction for w with the latest empirically determined 
corrections Ati? to the standard Newtonian/Einsteinian secular precessions of the perihelia 
of several solar system planets estimated by different teams of astronomers. We inferred 
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constraints on the minimum distance dx at which X may be located as a function of its 
echptic latitude /3x and longitude Ax for different values of its mass nix- We adopted the 
masses of Mars, 0.7 and 4 mjup. The use of the perihelia is effective in constraining 
dx for most of the sky locations: the availability of Aw for more than one planet allowed 
to mutually overlap the forbidden regions, thus further strengthening the lower bounds on 
rx also in those positions of the sky in which the use of a single planet at a time would, 
instead, yield weaker constraints. The most stringent bounds came from the extra-rates 
of the perihelia of Saturn, Mars, and, especially, the Earth thanks to the inclusion of the 
Jupiter VLBI data. Future improvements for Saturn, due to continuous ranging to Cassini, 
should yield tighter constraints. 

The minimum distance at which a remote body X with 70% of the mass of the Earth 
can exist is of the order of 250 — 450 au. For other values of the mass of X, we have 
150 — 200 au (Mars), and 3500 — 4500 (4 mj^p). Such lower bounds, which are independent 
of any speculation concerning the physical properties of X because they are solely based on 
its gravitational perturbations exerted on known planets, are tighter than those previously 
obtained by us in literature. Moreover, they are useful to better interpreting the results 
of past surveys aimed to directly detect X, and also to clarify what could realistically be 
expected from ongoing and future observational campaigns. 

We intend to further extend our analysis in future by using the newly determined 
co rrections to the stan dard secular precessions of the planetary perihelia and nodes as well 



bv lFienga et all fl201lf ). 



We remark that the corrections to the perihelion precessions used here were 
independently obtained by various astronomers without explicitly modeling the dynamical 
action of X itself. Thus, it might, in principle, have been somewhat "absorbed" and partially 
removed in the usual process of estimation of the solution's parameters like, e.g., the initial 
state vectors of the planets, especially if its magnitude was very small. To circumvent such 
an issue one should, actually, re-process the entire observational data set with modified 
dynamical force models explicitly including X itself, and solving for one or more dedicated 
parameters as well. Then, it would, e.g., be possible to look at the consequences of the 
inclusion of X on the values of the other usually estimated parameters. 
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